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We study the out-of-equilibrium dynamics of quantum systems with long-range interactions. Two different 
models describing, respectively, interacting lattice bosons and spins are considered. Our study relies on a 
combined approach based, on one hand, on accurate many-body numerical calculations and, on the other hand, 
on a quasi-particle microscopic theory. For sufficiently fast decaying long-range potentials, we find that the 
quantum speed limit set by the long-range Lieb-Robinson bounds is never attained and a purely ballistic behavior 
is found. For slowly decaying potentials, a radically different scenario is observed. In the bosonic case, a 
remarkable local spreading of correlations is still observed, despite the existence of infinitely fast traveling 
excitations in the system. This is in marked contrast with the spin case, where locality is broken. We finally 
provide a microscopic justification of the different regimes observed and of the origin of the protected locality 
in the bosonic model. 
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It is common wisdom that the propagation of a signal 
through a classical medium presents a distinct notion of 
causality, characterized by the progressive time growth of the 
spatial region explored by the signal. In spite of the intrin¬ 
sically non-local nature of the quantum theory, this familiar 
notion of locality is preserved in a wide class of quantum 
systems with short-range interactions. A milestone example 
is provided by the Lieb-Robinson (LR) bounds, which set a 
ballistic limit to the propagation of information, with expo¬ 
nentially small leaks outside the locality cone [1, 2]. The 
existence of LR bounds has many fundamental implications 
for thermalization, entanglement scaling laws, and informa¬ 
tion transfer in quantum systems [3]. A renewed interest in 
these topics is currently sparked by the impressive progress in 
the time-dependent control of ultracold-atom systems. Direct 
observation of cone spreading of correlations was reported in 
Refs. [4, 5]. 

The extension of the notion of locality to quantum systems 
with long-range interactions constitutes a fundamental chal¬ 
lenge. The paradigmatic model of long-range interactions 
considers an algebraic decay of some coupling term of the 
form V(R) ~ 1 /R a [6-10]. It applies either to the exchange 
coupling term in spin systems, as realized in cold ion crystals 
[11, 12], or to the two-body interactions in particle systems, 
as realized in ultracold gases of polar molecules [13], mag¬ 
netic atoms [14], and Rydberg atoms [15]. A remarkable fea¬ 
ture of long-range systems is that instantaneous propagation 
of information, in violation of locality, can take place when 
the exponent a is smaller than some threshold. This possi¬ 
bility is supported by the known extensions of the LR bounds 
to long-range interactions [16-18]. The latter yield “quasi- 
local” super-ballistic bounds for a> d, where d is the dimen¬ 
sion of the system, whereas for a < d no known generalized 
bounds exist, hence suggesting the breaking of quasi-locality. 
Evidence of the breaking of quasi-locality in one-dimensional 
(ID) Ising spin systems has been reported theoretically [6, 7] 
and experimentally in cold ion crystals [12, 19]. However, 
many questions remain open. For instance, although the ob¬ 


servations are compatible with the known long-range bounds, 
the propagation was found to be much slower than expected 
[7]. Hence, the bounds are usually not saturated and it is not 
clear that they provide a universal criterion for the breaking of 
quasi-locality. Moreover the threshold value for the breaking 
of locality in these systems is debated, and contrasting results 
have been put forward [6, 7]. To progress on these questions, 
it is of crucial importance to provide a unified understanding 
of a wider class of systems and, at the same time, to under¬ 
stand the microscopic origin of the breaking of quasi-locality. 

In this Rapid Communication, we study the out-of- 
equilibrium dynamics, induced by an interaction quench, of 
homogeneous ID quantum systems with long-range alge¬ 
braic interactions. We consider two different models, namely 
the long-range transverse Ising (LRTI) and Bose-Hubbard 
(LRBH) models. A quantitative analysis of quasi-locality 
in these systems is realized upon studying the earliest times 
at which information arrives at some fixed distance. On 
one hand, we perform ab-initio quantum many-body calcu¬ 
lations based on the time-dependent variational Monte Carlo 
(t-VMC) approach [20] . On the other hand, we provide a uni¬ 
fied analytical framework based on quasi-particle (QP) analy¬ 
sis. Both approaches consistently show that the two systems 
behave dramatically different. For a > 2 the LRTI model 
shows ballistic spreading of correlations with exponentially 
small leaks in time, hence leading to a strong form of quasi¬ 
locality. For 1 < a < 2 quasi-locality is still found. However, 
algebraic leaks in time appear, which can be traced back to a 
divergent group velocity at low momenta. For a < 1 quasi¬ 
locality is instead completely broken. This effect is traced 
back to the divergences of both the QP energy and velocity, 
which induce infinitely fast oscillations and a response time 
scale that goes to zero with the system size. Conversely, for 
the LRBH model, we find ballistic spreading of correlations 
for any value of a, analogous to what found in the short-range 
Bose-Hubbard model and in marked contrast with expecta¬ 
tions based on the lack of LR bounds for these systems. This 
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FIG. 1: Color online: Correlation spreading in long-range spin and boson models for various values of a. (a) Connected spin-spin correlation 
function in the LRTI model for the quench VJ- = h/2 —» Vf = /i/10. (b) Connected density-density correlation function in the LRBH model for 
the quench Ut = Vi=J^Uf = Vf = J/ 4 . Results obtained using the t-VMC approach for systems of L = 400 sites (for visibility, only a part 
is shown). The length unit is the lattice spacing and the time units are ft/h for (a) and h/J for (b). 


effect is traced back to the fact that the QP energy remains 
finite, which cancels non-local contributions for any value of 
a. All the observed regimes are explained by the unified QP 
analysis and shed new light on the microscopic origin of lo¬ 
cality in long-range quantum systems. 

Long-range transverse Ising model- We start with the 
long-range transverse Ising (LRTI) model, whose Hamilto¬ 
nian reads 

_ y _ erf erf 

^ + y ^ | •_ f ’ W 

i z #j\ l J\ 

where of, of are the Pauli matrices, h is the transverse field, 
V is the strength of the long-range spin exchange term, and in 
the following we set fi = 1 for convenience. Hamiltonian (1) 
is the prototype for long-range interacting quantum systems 
[6, 7, 21]. Moreover, it is experimentally implemented in cold 
ion crystals [22]. Evidence of the breaking of quasi-locality 
in information spreading has been reported for the ID LRTI 
model for sufficiently small exponents a [6, 7,12, 19], consis¬ 
tently with the absence of a long-range LR bound for a < 1. 
It was pointed out, however, that a model-dependent form of 
quasi locality may occur for specific initial states [7] with a 
complete understanding of the possible scenarios being de¬ 
bated. 

Asymptotically reliable results to reveal quasi-locality re¬ 
quire sufficiently long propagation times and sufficiently large 
systems. This is particularly crucial to determine precisely 
the nature of the dynamical regimes. To achieve this goal, we 
compute the unitary evolution of the correlation functions by 
means of the t-VMC approach [20, 23] (see Suppl. Mat.). The 
latter permits to simulate the dynamics of correlated quantum 
systems with an accuracy comparable to tensor-network meth¬ 
ods and proved numerically stable for unprecedented long 


times and large sizes. 

In the t-VMC calculations, we use a Jastrow wavefunction 
with long-range spin-spin correlations at arbitrary distance 
[24]. In order to avoid misleading finite size effects, which 
are usually strong in these issues [25], periodic boundary con¬ 
ditions (PBC) are used. For a lattice of size L with PBC, the 
interaction potential is taken as the sum of the contributions 
resulting from all the periodic images of the finite system. The 
Fourier components of the effective interaction potential are 

COS ( kfl ) 

then P{k) = 2 JT =1 --— = 2C1 a (k), where we have used 

n a 

the Poisson summation formula over the periodic images, k 
is an integer multiple of 2iz/L, and Cl a (k) is the Clausen co¬ 
sine function. In order to have a well-behaved potential in the 
thermodynamic limit, we set P (k = 0) = 0. It is the equiva¬ 
lent of the standard regularization procedure ensuring charge 
neutrality in the presence of electrostatic interactions [26]. 

We consider global quenches of the strength of the long- 
range interaction, Vi Vf. The results for the time- 
connected average G GG (R,t) = G GG (R,t) — G GG (R, 0) of the 
spin-spin correlation function G gg (R) = (ofof+a) are shown 
in Fig. 1 (a). We find three qualitatively different regimes. 
For a < 1, Fig. 1 (al), the propagation of correlations takes 
place on extremely short time scales and no cone-like struc¬ 
ture emerges. This is the signature of an efficient micro¬ 
scopic mechanism leading to the breaking of locality in the 
system. For a > 2, Fig. 1 (a3), a correlation cone with a 
well-determined velocity v clearly emerges. It is marked by 
a strong suppression of leaks in the region defined by R/t > v 
and space-time oscillations in the region R/t < v. In the in¬ 
termediate regime where 1 < a < 2, Fig. 1 (a2), a correla¬ 
tion cone-like structure is still visible but prominent non-local 
leaks are also appearing. 

In order to quantify more precisely the time-dependence 
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FIG. 2: (Color online) Panels (al,bl): Time-integrated spin-spin cor¬ 
relation functions G GG (R, t ) for the same data as in Fig. 1 (a). Super¬ 
imposed lines show the activation time t*(e) (see text), for £ ranging 
from 2 x 10 -2 (lighter lines) down to 5 x 10 -3 (darker lines). Panels 
(a2,b2): Behavior of the fitted exponent ) 3, computed within linear 
spin-wave theory, as a function of the cutoff parameter £ and for the 
system size L = 2 14 . The length unit is the lattice spacing and the 
time unit is h/h. 


of the leaks in the quasi-local regimes we study the 
time-integrated absolute value of the correlation function 
G° a (R,t) = \ /,' dt'\G GG (R.t')\. While it retains all the fea- 
tures of the signal propagation, it is less sensitive to time os¬ 
cillations. In Fig. 2 (al) we show the behavior of G° a (R,t) 
for a = 3. It clearly shows the sharp boundary of a ballis¬ 
tic cone. This is further assessed introducing a small cutoff 
£ and computing the first propagation time t*{R) such that 
G^ a (R,C) > £. The result is almost independent of £ and 
we find the scaling vt* = R P , with finite v and /3 ~ 1 to very 
good precision and up to large system sizes and long prop¬ 
agation times. The presence of a ballistic spreading, with 
exponentially suppressed leaks in time, is a stronger realiza¬ 
tion of locality than what expected from the looser long-range 
LR bound, which instead allows for polynomially suppressed 
leaks in time [30]. For 1 < a < 2 the same analysis of the 
leaks, shown in Fig. 2 (bl), reveals instead that polynomial 
leaks in time appear with an exponent /3 a, and a velocity 
v that vanishes with £. This is compatible with the long-range 
LR bound [16]. Remarkably, the regimes we find here for a 
global quench are the same qualitative regimes that have been 
identified for a local quench in the LRTI model in Ref. [6]. 

Long-range Bose-Hubbard model- We now turn to the 
long-range Bose-Hubbard (LRBH) model, which describes 
interacting spinless bosons in a periodic potential with 
nearest-neighbor tunneling and long-range two-body interac¬ 
tions. The Hamiltonian reads 


(iJ) 


]bj + h.c)j + y («. " i) + ^ £ J 7 


¥J 
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where bi(bj) destroys (creates) a boson on site i, ni = b[bi 
is the particle number operator, J is the tunneling amplitude, 


U is the on-site interaction energy, V is the strength of the 
interaction potential, and we set h = 1 for convenience. The 
short-range (V = 0) case is now routinely realized in ultracold- 
atom experiments and the long-range (V ^ 0) case applies to 
polar molecules [13], magnetic [14] and Rydberg atoms [15]. 

We perform t-VMC calculations using a Jastrow wave- 
function incorporating density-density correlations at arbi¬ 
trary large distances [27] (see Suppl. Mat.). For simplic¬ 
ity we choose U = V, fix the density at half filling (n = ^), 
and consider the connected density-density correlation func¬ 
tion G™ (R) = ( nini+R) —n 2 . We study global quenches in the 
interaction strength Vi —>Vf. The results for various values 
of the exponent a are shown in Fig. 1 (b). Surprisingly, we 
find here that the LRBH model exhibits the same qualitative 
behavior for all values of a, in marked contrast with the LRTI 
model. Within numerical precision, we always find a purely 
ballistic cone spreading of correlations at some velocity v. 
The long-range LR bound is therefore never saturated and, for 
every value of the exponent a, the spreading is qualitatively 
identical to the short-range case. Hence, in the LRBH model 
quasi-locality appears to be strongly protected even for very 
long-range interactions. This is further confirmed by a pre¬ 
cise analysis of the leaks, along the same lines as for the LRTI 
model. It always yields a scaling of the form vt* = R P with 
/3 = 1 and the signal is exponentially suppressed for times out 
of the locality cone, i.e. when t < R/v. 

Quasi-particle analysis- The radically different behav¬ 
iors of the LRTI and LRBH models are particularly striking 
because they share the same class of long-range interactions 
and are therefore subjected to the same universal long-range 
LR bounds [16] . In order to understand the different behaviors 
of the two models, at a miscroscopic level, we use a general 
QP approach. The latter has a broad range of applications, e.g. 
universal conformal theories [28], spin systems [6], superflu¬ 
ids [29], and Mott insulators [4]. A generic time-dependent 
two-body, connected correlation function in a translation in¬ 
variant model with well-defined QP excitations can be written 
as 


/ K dk ( 

—&(k) {cos m+ 


where E^ is the ^-momentum QP energy of the post-quench 
Hamiltonian and & (k) is the weight associated to each QP. 
This general form states that the collective excitations of the 
system are coherent superpositions of pairs carrying excita¬ 
tions of momentum k and traveling in opposite directions. 
Whereas e[ depends only on the post-quench Hamiltonian, 
in general & (k) instead depends both on the pre- and post- 
quench Hamiltonians [31]. 

In the LRTI model and in the regime of large transverse 
field (h V), considered in the t-VMC calculations, we can 
apply linear spin wave theory [11]. The QP energy and weight 
read, respectively, e£ = 2 ^/h[h-\-VP(k)] and & GO (k) = 


cos (^kR — 2E^t^ + cos (kR + 2 E^t^ | 
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FIG. 3: (Color online) QP analysis of the group velocity for the 
LRBH model with exponent a = 1/2, density n = 1/2 and U = V = 
J / 4. In the inset, the stationary phase weights corresponding to the 
three solution branches of the equation 2v g (k) = v are shown. Wave- 
vectors are in units of the inverse lattice spacing and velocities are in 
units of the hopping amplitude J. 


2P(k)(Vi-V f ) 

E[[\+P{k)V f /h\ 


, where P(k) are the Fourier components 


of the interaction potential (see Suppl. Mat.). Let us then an¬ 
alyze the outcome of Eq. (3). 

The ballistic behavior observed for a > 2 can be understood 
from stationary phase analysis. Along the line R = vt, it yields 
the dominant contribution 


G^(RJ) 


&(k*) 



t 


cos (k*R - 2 e[j + </>) , (4) 


where k* is solution of v = 2v g (k*), v g is the group velocity, 
and 0 is a time-independent phase. For a > 2, the group ve¬ 
locity is bounded and we find a ballistic cone spreading of cor¬ 
relations with a velocity that is given by twice the maximum 
group velocity. We have checked that it agrees well with the 
numerics where a cone propagation is also found. A quantita¬ 
tive analysis of the leaks within the spin-wave approach con¬ 
firms that they are exponentially suppressed in time. More¬ 
over, the exponent /3 is found to smoothly approach unity 
upon reducing the cutoff parameter £ [see Fig. 2 (a2)], in 
agreement with the t-VMC results. 

For 1 < a < 2 the group velocity is instead unbounded, 
since it exhibits the infrared divergence v g (k) °c k~ l 2-a L 
Hence the correlation front is no longer given by a well- 
determined velocity but by the coherent superposition of in¬ 
finitely fast modes at low momenta. More precisely, the be¬ 
havior of the leaks can be found from the asymptotic R °o 
expansion of Eq. (3), retaining only the contributions of the 
divergent velocities. For a = 3/2, it can be computed ex¬ 
actly, and yields G GG {R °°,t) ~ F(t) x t/R a , where F(t) is 
a bounded function of time. This scaling agrees with and ex¬ 
plains the t-VMC results. It is further confirmed by the analy¬ 
sis of the leaks within the spin-wave approach. The exponent 


jS is found to smoothly approach the interaction exponent a 
upon reducing the cutoff parameter £ [see Fig. 2 (b2)]. The 
same was found for other values of a between 1 and 2. 

For a < 1, both the QP energy and the group velocity di¬ 
verge for k 0, respectively, as e[ = eok~ I “2“ I, v g (k) 

k ~I “2“ I and eo = 2yJhfVf. In particular, it is the energy di¬ 
vergence which sets the breaking of quasi-locality in this case. 
The latter gives rise to a rapidly oscillating factor of the form 

cos (eot/k^ ^ in Eq. (3), which leads to a non-analytic t m 0 
delta-kick in the thermodynamic limit. More precisely, an ex¬ 
act asymptotic expansion of the correlation function (3) can 
be derived in the limit of small propagation times t and large 
distances R. Keeping the relevant small quasi-momenta, it 
yields 


G GG (R,t-> 0)~ limA- 

L—>oo 


sin 


(i“«) 


eo t 


cos (R/L) +B (eot) 


R 2-a 


R 


1+q ’ 


( 5 ) 

where A and B are finite numerical constants, which depend 
on the microscopic parameters of the model. A remarkable 
consequence of this expression is that the first term yields 
an instantaneous contribution to the signal, on a time scale 


(1-q) 

eoT =l/L 2 , independent of the distance R and with an ex¬ 
ponent set by the divergence of the QP energy. This implies 
that the system reacts on a time scale inversely proportional 
to the system size, yielding an efficient mechanism for the 
breaking of locality. 

An analogous microscopic analysis can be carried out for 
the LRBH model in the weakly interacting superfluid regime, 
considered in the t-VMC calculations. In this regime the 
quantities e[ and ^ (k) are found by Bogoliubov analysis 

[29], which yields E^ = J £^ [e* + 2n (Uf -\-V/P (k))\ and 


^ ’ where £k = 4/81112 
is the free-particle lattice dispersion and n is the particle den¬ 
sity (see Suppl. Mat.). 

For a > 1, the origin of the observed ballistic behavior is 
traced back to the fact that QP velocities are bounded. The 
propagation of correlations is dominated by the stationary- 
phase points of Eq. (4) and the correlation cone velocity is 
given by twice the maximum group velocity. 

For a < 1, the group velocity diverges as v g (k —» 0) oc 

I 1—q I 

k \ pT I, whereas the QP energy is always finite. Hence, at 
variance with the LRTI model, the correlation function does 
not exhibit any instantaneous kick at t = 0 such as that of Eq. 
(5) and quasi-locality is preserved. Moreover, although this 
case is formally analogous to the intermediate regime, with 
polynomial leaks in time, found for the LRTI model, a purely 
ballistic spreading is instead found within numerical precision 
in the LRBH model. To understand this, let us come back to 
the stationary-phase approach of Eq. (4). Due to the specific 
form of the group-velocity dispersion in the LRBH model, 
shown in Fig. 3, the equation v = 2v g (k k ) has up to three sep¬ 
arate solutions for a given velocity v. The correlation function 
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is thus dominated by three contributions (I, II, and III ) of the 
form of Eq. (4). The behavior of the corresponding weights, 

W (v) = & nn (k*) / yj\d%Ef (k*)/j\, along the three branches 
is shown in the inset of Fig. 3. The largest weights, corre¬ 
sponding to the velocities dominating the propagation, belong 
to the regular branches (II and III). The latter extend up to 
a certain maximum velocity v max , which effectively sets the 
correlation cone velocity v c . The infinitely fast modes, corre¬ 
sponding instead to branch /, have a weight which is polyno¬ 
mial^ suppressed at large velocities. The exponent for the al¬ 
gebraic decay can be derived from the known small k behavior 

of P(k), which leads to W(v °o) oc v 1 2 ( 1 -«) I. For v ~ v max , 
the weights of these modes are already few orders of magni¬ 
tude smaller than the quasi-local, ballistic modes. This separa¬ 
tion of scales is responsible of the effective suppression of the 
infinitely fast non-local modes in the FRBH model. Notice 
that in the FRTI model the irregular branch of the infinitely 
fast modes is never protected by a finite-velocity branch. This 
is a consequence of the monotonic behavior of the QP group 
velocity in the spin model, yielding a single branch of solu¬ 
tions for the stationary phase equation. 


Discussion- In summary, the radically different behav¬ 
iors of the FRTI and FRBH models, found both in t-VMC 
and QP analysis, show that specific microscopic properties of 
the system, not accounted for in universal FR-like bounds, 
play a major role in quasi-locality. This result sheds new light 
on the dynamics of long-range quantum systems. Yet many 
questions remain open and are worth being investigated in the 
future. For instance, it is expected, on the basis of univer¬ 
sal bounds, that the critical exponents for the breaking of lo¬ 
cality depend on the system dimension [16]. It would thus 
be of utmost interest to study the counterparts of the effects 
discussed here in dimension higher than one, which could 
be done by a straightforward application of the present ap¬ 
proach [23]. Moreover, it would be interesting to study the 
same FRTI and FRBH models in a regime of stronger interac¬ 
tions, where an extension of the t-VMC and QP analysis tak¬ 
ing into account relevant excitations can be developed. In this 
regime, the FRBH model maps onto an effective spin model 
and might therefore exhibit radically different properties than 
those found in this work. 


We acknowledge discussions with M. Cheneau, M. Fagotti, 
M. Holtzmann, D. Porras, and F. Tagliacozzo. This re¬ 
search was supported by the European Research Council 
(FP7/2007-2013 Grant Agreement No. 256294), Marie Curie 
IEF (FP7/2007- 2013 - Grant Agreement No. 327143), 
FET-Proactive QUIC (H2020 grant No. 641122). It was 
performed using HPC resources from GENCI-CCRT/CINES 
(Grant c2015056853). Use of the computing facility cluster 
GMPCS of the FUMAT federation (FR FUMAT 2764) is also 
acknowledged. 


[1] E. H. Lieb and D. W. Robinson, Comm. Math. Phys. 28, 251 
(1972). 

[2] J. Eisert, M. Friesdorf, and C. Gogolin, Nat. Phys. 11 , 124 
(2015). 

[3] B. Nachtergaele and R. Sims, IAMP News Bulletin , 22 (2010). 

[4] M. Cheneau, P. Barmettler, D. Poletti, M. Endres, P. Schausz, 
T. Fukuhara, C. Gross, I. Bloch, C. Kollath, and S. Ruhr, Nature 
481 , 484 (2012). 

[5] T. Langen, R. Geiger, M. Kuhnert, B. Rauer, and J. Schmied- 
mayer, Nat. Phys. 9, 640 (2013). 

[6] P. Hauke and L. Tagliacozzo, Phys. Rev. Lett. Ill, 207202 
(2013). 

[7] J. Eisert, M. van den Worm, S. R. Manmana, and M. Kastner, 
Phys. Rev. Lett. Ill, 260401 (2013). 

[8] D. Vodola, L. Lepori, E. Ercolessi, A. V. Gorshkov, and 
G. Pupillo, Phys. Rev. Lett. 113 , 156402 (2014). 

[9] J. Schachenmayer, A. Pikovski, and A. M. Rey, New J. Phys., 
065009 (2015). 

[10] M. A. Baranov, M. Dalmonte, G. Pupillo, and P. Zoller, Chem. 
Rev. 112 , 5012 (2012). 

[11] X.-L. Deng, D. Porras, and J. I. Cirac, Phys. Rev. A 72, 063407 
(2005). 

[12] P. Richerme, Z.-X. Gong, A. Lee, C. Senko, J. Smith, M. Foss- 
Feig, S. Michalakis, A. V. Gorshkov, and C. Monroe, Nature 
511 , 198 (2014). 

[13] H.-P. Buchler, G. Pupillo, A. Micheli, and P. Zoller, “Con¬ 
densed matter physics with cold polar molecules,” in Cold 
Molecules (2009). 

[14] C. Menotti, M. Lewenstein, T. Lahaye, and T. Pfau, in AIP 
Conference Proceedings , Vol. 970 (AIP Publishing, 2008) p. 
332. 

[15] T. Macri and T. Pohl, Phys. Rev. A 89, 011402 (2014). 

[16] M. B. Hastings, (2010), arXiv: 1008.5137. 

[17] Z.-X. Gong, M. Foss-Feig, S. Michalakis, and A. V. Gorshkov, 
Phys. Rev. Lett. 113 , 030602 (2014). 

[18] M. Foss-Feig, Z.-X. Gong, C. W. Clark, and A. V. Gorshkov, 
Phys. Rev. Lett. 114 , 157201 (2015). 

[19] P. Jurcevic, B. P. Lanyon, P. Hauke, C. Hempel, P. Zoller, 
R. Blatt, and C. F. Roos, Nature 511 , 202 (2014). 

[20] G. Carleo, F. Becca, M. Schiro, and M. Fabrizio, Sci. Rep. 2, 
243 (2012). 

[21] M. van den Worm, B. C. Sawyer, J. J. Bollinger, and M. Kast¬ 
ner, New J. Phys. 15 , 083007 (2013). 

[22] D. Porras and J. I. Cirac, Phys. Rev. Lett. 92, 207901 (2004). 

[23] G. Carleo, F. Becca, L. Sanchez-Palencia, S. Sorella, and 
M. Fabrizio, Phys. Rev. A 89 , 031602 (2014). 

[24] F. Franjic and S. Sorella, Progr. Theor. Phys. 97, 399 (1997). 

[25] J. Haegeman, C. Lubich, I. Oseledets, B. Vandereycken, and 
F. Verstraete, (2014), arXiv: 1408.5056. 

[26] N. W. Ashcroft, Solid State Physics (Holt Rinehart & Winston, 
1976). 

[27] M. Capello, F. Becca, M. Fabrizio, and S. Sorella, Phys. Rev. 
Lett. 99 , 056402 (2007). 

[28] P. Calabrese and J. Cardy, Phys. Rev. Lett. 96 , 136801 (2006). 

[29] S. S. Natu and E. J. Mueller, Phys. Rev. A 87, 053607 (2013). 

[30] We focus here exclusively on the time-dependence of the leaks, 
and we find regimes with exponentially suppressed leaks in 
time. Notice that the spatial dependence of the correlation func¬ 
tion, even for short-range Hamiltonians, can instead exhibit al¬ 
gebraic behavior, associated to quasi-long-range order in ID. 

[31] The validity of the QP picture in this context is corroborated 





6 


by the good quantitative agreement between the QP dynamics 
given by Eq. (3) and the t-VMC results (see Suppl. Mat.). 


